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I .  Introduction 

Coherent  transient  interactions  between  materials  and  optical  fields  are 

one  of  the  most  interesting  subjects  in  physics.  Investigations  of  transient 

optical  effects  have  become  very  important  in  the  development  of  coherent 

optical  spectroscopy  and  fast  response  nonlinear  optical  materials.  Such 

studies  cover  such  a  wide  spectrum  of  subjects  and  materials  as  photon  echo, 

free  induction  decay  and  optical  nutation  in  gaseous  materials,^  ultrafast 

dynamical  blue  shift  of  the  exciton  resonance  accompanied  by  strong  bleaching 

6-9 

in  quantum  wells  in  heterostructure  semiconductors,  nonlinear  optical 

response  of  excitons^  ’  ^  and  large  susceptibilities^’^  in  polymers,  optical 

nutation  in  direct-gap  bulk  semiconductors , ^ ^  and  so  forth. 

A  nonlinear  optical  material  can  be  defined  as  one  whose  optical 

properties  change  when  light  shines  upon  it.  In  recent  years,  the  study  of 

coherent  nonlinear  optical  processes  produced  by  the  excitation  of 

6-9 

semiconductors  in  the  transparent  region  well  below  the  absorption  edge  has 

been  extended  to  polymers . ^ ^  When  the  frequency  of  the  laser  beam  falling 

on  semiconductors  is  tuned  below  the  exciton  resonance,  it  is  known  that 

virtual  excitons  are  generated.  These  excitons  interact  with  photons  in 

7  9-11 

exactly  the  same  way  as  real  excitons.  ’  Such  excitonic  effects  can  be 

9 

explored  more  naturally  by  probing  the  susceptibility  and  looking  for 
transient  changes  that  persist  as  long  as  the  excitation  ^e  imaginary  part 
of  the  susceptibility  is  usually  examined  via  the  absorptic  .  oectrum,  and  its 
real  part  is  explored  by  measuring  the  index  of  refraction,  provided  that  the 
sample  is  strongly  excited  well  below  the  conduction  band  edge. 

It  is  well  known  that  polymers  exhibit  large  nonlinear  optical 
susceptibilities.  In  comparison  with  other  materials,  polymers  have  the  added 
advantage  of  having  extremely  fast  ground-state  recovery  times,  which  is  why 


they  attract  much  attention  for  their  potential  applications  in  future 

generations  of  signal  processing  devices.  Polydiacetylene  (PDA)  may  be  a  good 

(  3  ) 

candidate  for  such  considerations  because  of  its  large  x  and  small 

transmission  loss  a  which  can  be  reduced  to  as  low  as^  a  <  1  db/cm.  Thus  the 

(3) 

material  possesses  a  fairly  large  ratio  x  /a<  which  is  usually  the  measure 
of  the  usefulness  of  a  material  employed  in  switching  devices.  Besides,  PDA 
is  easier  to  use  in  the  construction  of  waveguides . ^  In  spite  of  the 
considerable  amount  of  research  carried  out  so  far,  however,  the  mechanism 
that  governs  its  nonlinear  optical  response  is  still  not  well  understood, 
although  tremendous  progress  has  made  it  ripe  for  further  examination. 

In  this  paper,  we  study  the  transient  behavior  of  the  induced  optical 
susceptibility  of  polymers  due  to  irradiation  by  an  ultrashort  pulse  of  a  pump 
laser  field.  We  take,  as  an  example,  polydiacetylene- toluene  sulfanate  (PTS) 
which  is  irradiated  by  a  strong  pump  laser  and  a  weak  probing  laser.  The  two 
lasers  may  have  pulse  duration  of  femtoseconds  to  picoseconds.  For  the  PTS 
system  we  adopt  the  two -level  model^’  to  describe  its  ground- state  and 

virtual  exciton  state.  The  lifetime  of  the  exciton  and  the  mediating  effects 
of  the  phonon  are  introduced  phenomenologically.  We  use  the  perturbation 
technique  to  find  the  differential  nonlinear  susceptibility  by  solving 
analytically  to  all  orders  in  the  pump  field  and  to  first  order  in  the  test 
field. 

II.  Two-level  model 

The  two  electronic  states  of  PTS  are  denoted  by  |+>  and  |->.  We  assume 
a  dipole  transition  and  define  the  dipole  operators 


S'  -  |  -x+ 


(la) 


S+  -  j+X- 


(lb) 


with  the  inversion  operator  given  by 

SZ  =*  2  (  l+X+l  -  |  -X-  |  )  .  (lc) 

For  a  single  virtual  exciton  and  a  single  phonon  mode,  we  can  write  the  free 
Hamiltonian  as 


H  =■  Yw  S2  +  tfw.btb.  ,  (2a) 

O  X  111 

where  bt  (b.)  creates  (annihilates)  a  phonon  in  mode  i  with  frequency  cj  and 
i  l  1 

11 

w  is  the  exciton  frequency.  For  the  interaction  Hamiltonian  we  consider 
that  each  virtual  exciton  interacts  with  both  a  strong  pump  field  of  amplitude 
Ep  and  frequency  and  a  weak  test  field  of  amplitude  Et  and  frequency  . 
Since  the  test  field  is  very  weak,  we  assume  the  mediating  effects  of  the 
phonons  on  the  dipole  transition  to  be  induced  only  by  the  pump  field,  as 
shown  in  Fig.  1.  Therefore,  we  can  write  the  interaction  Hamiltonian  in  both 
the  rotating-wave  approximation  and  the  interacticn  representation  as 


,  lVt  , 

HP1  -  - 1  Av  1  bIs' +  h  c- 


(2b) 


for  the  interaction  between  the  pump  field  and  the  PTS ,  and 


1 

"c  -  -  2  "Ete 


S  +  h.c. 


(2c) 


for  the  interaction  between  the  test  field  and  the  PTS ,  where  A  =-  w  +  w 

Pj_  P  i 

^x>  -  u>x,  n  is  the  dipole  matrix  element,  and  A  is  a  coupling  constant 

proportional  to  the  transition  dipole  moment  p.  Here  we  have  neglected  the 

momentum  dependence  of  virtual  excitons^ '  ^  and  the  damping  effect  of  the 

phonons.  This  Hamiltonian  uses  dipole  or  projection  operators  to  describe  the 

dipole  transition  between  the  ground  state  and  the  excitonic  state  of  the 

material  by  external  fields.  It  is  different  from  that  of  Ref.  11,  where 

creation  and  annihilation  operators  were  used  for  virtual  excitons. 

For  simplicity,  we  assume  that  each  phonon  mode  has  only  two  states,  | 1> 

and  |0>,  and  that  the  two  and  higher-phoncn  eff.ects  are  negligible.  The 

initial  condition  we  consider  here  is  that  the  system  electronically  is  in  its 

ground  state  with  no  virtual  exciton  present,  while  the  phonon  mode  is  in  its 

excited  state  |l>.  We  define  Ci+(t)  as  the  probability  amplitude  of  the  state 

|±>  when  there  is  the  i-th  phonon  mode  only,  of  which  the  time  derivative  C^+ 

is  composed  of  three  parts:  pump  field- induced  part  Cf+,  test  field- induced 
•  t  •  d 

part  Ch  + ,  and  damping  part  C^+.  By  using  the  Schrodinger  equation  and  the 

14  15 

interaction  Hamiltonians  (2),  we  obtain  ’ 


6 


(t) 


'7Ci+ 


(4a) 


•  d 

c“_(t)  -  0 


(4b) 


where  7  is  the  damping  rate  describing 

•  •  •  P 

exciton.  Therefore,  from  C.,  -  C .  . 

i±  i± 

equations  for  Ch  +  : 


the  short  lifetime  of  the  virtual 

•  t  •  d 

+  C. ,  +  C.,  we  have  the  following 
i±  i±  & 


c. 

1+ 


(t) 


-  iA 


iQ*  e 
P 


in*  e 

t 


-  iA  t 
t 


C. 

l- 


7C, 


1+ 


(5a) 


C._(t) 


C. 

1+ 


+  in 


iAtC 


Ci+ 


(5b) 


Since  the  test  field  is  in  most  cases  much  weaker  than  the  pump  field, 
we  can  safely  apply  perturbtaion  theory  to  solve  (5) .  The  zeroth-order 
equations  are  obtained  by  setting  »  0  in  (5),  namely, 


in*  e 
P 


-  iA  t 
Pi 


(6a) 


(6b) 


In  what  follows,  we  assume  for  simplicity  n  -  p* ,  A 
With  the  initial  conditions  C°+(0)  -  0  and  C°  (0)  -  1, 


-  A*  and  n  -  0* 

p,  t  p,t 

we  find  the  solutions 


sinfl.t  exp[-(iA 
1  P4 


+  7)t/2] 


(7a) 


7  -  iA  (iA  -7) c/2 

p .  p . 

o  r  1  1 

CL  -  (cosOt  +  — — -  sinO^t)  e 

i 


where  we  have  defined  the  complex  Rabi  frequency 


1  2  2  4 

°i  "  2  [(AP.  +  +4V 


The  real  part  of  the  Rabi  frequency  is  given  by 


nR  =  [ (7x2  +  y2  +  X)/S]h 


which  determines  the  oscillation  frequency,  and  the  imaginary  part  by 


Oj  -  -[(7x2  +  y2  -  x)/8]4 


which  determines  the  envelope  of  the  oscillation.  In  (9)  we  have  defined 
2  2  2 

x  -  A  +40  -  7  and  y  -  2A  7,  with  A  <0. 

p.  p  J  p.  p. 

1  1  1 

The  solutions  up  to  first-order  perturbation  can  be  obtained  by 
inserting  the  zeroth-order  solutions  (7)  into  the  right-hand  side  of  (5) .  The 
results  are 


Ci±(t)  "  Ci±(t)  +  6Ci±(t)  +  Ai± 


exp[iA  t  -  (li  +  7)C] 

cr  _  - -i - 

6Ui-  0  1  2  2 

1  “i  [f(lAp  +  7)  -  i\]  +  0 ‘ 
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x 


[  t^CiA 


Pi 


+  7> 


iA  )  sinfl.t  +  fl.cosfl.t] 

t  i  li 


(11a) 


SC. 

i+ 


m  exPf-iV  +  2(iAp  -  -r)t] 

_ t  _  _  ri 


[iAc  -  2 ( i^p  *  7)1 2  +  n[ 


2  9 

x  ([ (A_  +  i7)  +  20‘  +  iA  (iA  -  7)]  sinft.t 

Pf  p  t'  p  1,1  1 


2fli[i(Ac  -  Ap  )  +  7]  cosfl^t) 


(lib) 


The  integration  constants  A^  +  in  (10)  are  determined  by  the  initial 
conditions . 

The  expectation  value  of  the  dipole  moment  of  a  virtual  exciton  is 


■iu  t  -in  t 

u.  (t)  -  /i< s‘(t)>  e  X  -  nC*  C.  e  x 
1  ^  1-1+ 


o*  -iu  t 

/i(C  +  $Ct  +  A*  )(C?,  +  SC.  +  A.  )  e  X 

1  1 '  1.  *  It  It 


(12) 


It  is  evident  that  the  dipole  moment  induced  by  the  test  beam  is  given  by 


-  ^(Ci_*iCi+  +  C°+5C*_)e 


'  iu  t 
x 


(13) 


Therefore,  the  nonlinear  optical  susceptibility  experienced  by  the  test  beam 
is  simply 


9 


-lw  t 

XL  -  ni5/ii(t)/(Ete  C  ) 


2  Ap  '  i7 

in  e‘7t  (— 7  [cosOtt  +  i  — ^  sinfltt] 

L  2*  20.  z,1  1  2ni 

l  1 


2  2 

x  [ ( (A  +  i7)  +20  -  A  (A  +  i7) }  sinO.t 

1  P  t  p.  i 


20. ( i(A  -  A  )  +  7>  cosO.tj 
L  C  pi  1 


0 


*77  »ln01t((i(14 


•  7>  - 


l°ll  Z2 


iAc)  sinOtt  +  0*  cosOtt] }  ,  (14) 


where 


'H, 


7) 


2  2 

iAj  +  0. 


(15a) 


t2(iV 

l 


-  7) 


2  2 
iA  ]  +  0* 

t  l 


(15b) 


We  have  introduced  in  (14)  the  optically- induced  transient  virtual 

exciton  density  n.  which  is  taken  to  be  proportional  to  the  probability  of 

18 

virtual  exciton  state, 


n,  - 


n  C?  C?* 
s  1+  1+ 


0 


n  — ®  7  IsinO.tl2  exp(-yt) 

5  |ntl2 


(16) 


10 


where  we  have  assumed  that  all  the  virtual  excitons  are  induced  by  the  pump 

field  E  ,  since  those  induced  bv  the  weaker  probe  field  are  much  less  likely, 

P 

and  n^  is  the  saturation  density  of  the  virtual  excitons.  In  this  paper  we 

19 

only  consider  the  two  phonon  modes  coupling  most  strongly  to  the  exciton,  so 
that  we  have  the  nonlinear  susceptibility 

2 

■  <17) 

i-1 

The  transient  behavior  of  the  nonlinear  optical  susceptibility  is 

calculated  numerically  as  a  function  of  time  and  the  frequencies  u>  and  u>  . 
The  results  are  presented  and  discussed  in  the  next  section. 

III.  Results  and  discussion 

In  our  numerical  study,  we  have  used  the  unit  jrf  —  1 .  Other  parameters 

are  chosen  as^’^'^  -  0.184  eV,  “  0-258  eV,  7  -  0.02  eV,  ~|m|  “ 

77-  eV,  0  -  5  x  lo'4  eV  and  w  -  1.978  eV.  The  real  and  imaginary  parts  of 

40jt  p  x 

are  computed  separately.  Figure  2  shows  snapshots  at  different  instants  of 

time  of  the  variation  of  with  the  pump  detuning.  The  variation  with  time 

for  fixed  detuning  is  shown  in  Fig.  3.  Figure  4  shows  snapshots  at  different 

instants  of  time  of  the  imaginary  part  of  x  when  the  resonance  condition  is 

satisfied,  i.e.,  w  -  w  +  . 

x  p  2 

It  is  clearly  seen  from  Fig.  2  that  both  the  differential  absorption  and 
refraction  are  direct  effects  of  the  virtual  excitons  induced  by  the  pump 
field.  The  mediating  role  played  by  the  phonon  is  substantial  when  the 
duration  of  the  pump  pulse  is  long  enough  (see  the  two  peaks  of  the  imaginary 
part,  where  w  —  w  +  w.).  This  is  in  agreement  with  the  experimental  results 


11 


of  Ref.  11.  However,  when  the  duration  of  the  pump  pulse  is  short,  e.g.,  less 
than  the  lifetime  of  the  virtual  exciton,  the  peak  structure  of  the  imaginary 
part  is  not  clear,  i.e.,  we  cannot  see  evident  optical  phonon  -  mediated 
effects.  It  should  be  noted  at  this  point  that  the  present  results  are 
different  from  those  reported  in  Ref.  21  in  which  steady-state 
susceptibilities  are  calculated  in  the  three-level  two-field  model  without 
phonon  effects. 

Optical  nutation  can  be  observed  in  polymers  in  the  ultrashort  time 
regime  when  virtual  excitons  are  generated  by  a  pump  field  tuned  below  the 
exciton  resonance.  This  situation  can  be  seen  from  Fig.  3.  However,  when  the 
sum  of  the  pump  field  frequency  and  one  of  the  phonon  frequencies  is  on 
resonance  with  the  exciton  frequency,  the  Rabi  oscillation  amplitude  becomes 
smaller.  This  reflects  the  fact  that  the  optical  phonon  plays  a  crucial  role 
in  the  optical  nutation  of  PTS .  It  is  also  clearly  shown  that  sometimes  the 
oscillation  centers  of  Imxt  change  with  time.  This  is  quite  different  from 
nutation  phenomena  predicted  in  Refs.  14  and  15  for  semiconductors  where  the 
oscillation  amplitude  decreases  monotonically  like  damped  oscillators. 

From  Fig.  4,  we  can  see  the  transient  behavior  of  the  spectral  hole 

burning,  which  also  needs  sufficient  long  duration  of  the  pump  pulse  to  have 

an  effect.  The  hole  is  close  to  the  exciton  resonance,  and  there  is  a  deeper 

hole  when  w  -  w  +  u. .  This  indicates  that  there  is  energy  transfer  from  the 
x  p  1 

pump  field  to  the  test  field  via  optical  phonons  and  virtual  excitons,  also  in 
qualitative  agreement  with  experiments . ^ ^ 

Since  it  takes  time  for  the  pump  field  to  build  up  enough  virtual 
excitons  and  since  the  virtual  excitons  are  short-lived,  their  effects  are 
appreciable  and  stable  raos  -  T  the  time  during  the  pulse  (on  the  order  of 
picoseconds).  The  respon  vae  can  be  anywhere  between  20  -  200  femtoseconds 


according  to  our  calculation.  Both  the  absorption  and  refraction  parts  show 
up  almost  right  away  and  taper  off  rather  slowly  after  reaching  their  peaks. 
So,  for  longer  excitation  time  of  the  pump  field,  say,  several  times  of  the 
exciton  lifetime,  the  steady-state  method  may  be  reliable. 

In  addition,  we  have  also  found  that  increasing  the  intensity  of  the 

pump  field  leads  to  a  deeper  and  wider  hole  in  the  absorption  spectrum,  which 

is  qualitatively  in  agreement  with  the  steady-state  results  of  Ref.  11.  What 

seems  to  be  more  interesting  is  that  we  find  a  blue  shift  of  0.005  eV  in  our 

calculation.  This  represents  joint  effects  of  both  the  phonon- induced  optical 

Stark  shift  and  the  electronic  state  optical  Stark  shift.  Such  effects  have 

20 

been  observed  in  a  recent  experiment  in  which  the  differential  transmission 
of  polydiacetylene  -  toluene  sulfanate  has  been  measured  by  means  of  coherent 
inverse  Raman  spectroscopy.  It  is  demonstrated  in  Ref.  20  that  this  ac  Stark 
effect  is  important  in  the  determination  of  the  nonlinear  optical  response  of 
PTS  even  in  the  small  signal  limit.  The  effect  is  reproduced  theoretically  in 
Ref.  20  by  including  an  extra  term  to  modify  the  exciton  resonance  frequency 
in  the  simple  model  of  Schmitt-Rink  et  al.^  However,  we  need  not  do  so,  and 
the  treatment  in  both  Refs.  11  and  20  assumes  a  steady-state  solution  while 
the  experimental  conditions  are  more  transient- like . 
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Figure  captions 

1.  Two-level  model  representing  the  excitonic  state  | +>  and  ground  state  |->. 

The  virtual  excitons  are  induced  by  a  strong  pump  field  with  frequency  w  . 

Each  transitions  is  mediated  by  a  single  phonon  of  frequency  ,  such  that 

u  +  w.  —  w  on  resonance, 
p  1  x 

2.  Snapshots  of  the  real  (curve  1)  and  imaginary  (c  rve  2)  parts  of  the 

nonlinear  susceptibility  versus  pump  detuning  at  various  times  for  w  = 
wx:  (a)  t  -  25  (eV)"1,  (b)  t  -  100  (eV)'1,  (c)  t  -  500  (eV)"1,  (d)  t  - 

3500  (eV)'1. 

3.  Real  (curve  1)  and  imaginary  (curve  2)  parts  of  versus  time  for  =* 

w  :  (a)  w  -  1.6  eV;  (b)  w  -  1.72  eV. 

x  p  p 

4.  Snapshots  of  the  imaginary  part  of  the  susceptibility  versus  probe 
frequency  when  u  -  1.72  eV:  (a)  t-25  (eV)  \  (b)  t-  500  (eV)  ,  (c)  t 


-  3500  (eV) 


Ka.  1 
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